from sympy import *
# 输入系数矩阵A
A = Matrix([[1, 1, -3,-1], [3, -1, -3, 4], [1, 5, -9, -8]])
# 输入常数项矩阵b
b = Matrix([1, 4, 0])
# 利用逆矩阵求解方程组的一个特解
x_t = A.pinv()*b
# 计算线性方程组的基础解系
B = A.nullspace()
# 输出结果
print('方程组的一个特解为:\n',x_t)
print('方程组的基础解系为:\n',B)